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Abstract 

General properties of perturbed conformal field theory interacting with quantized Liou- 
ville gravity are considered in the simplest case of spherical topology. We discuss both short 
distance and large distance asymptotic of the partition function. The crossover region is 
studied numerically for a simple example of the perturbed Yang-Lee model, complemented 
in general with arbitrary conformal "spectator" matter. The latter is not perturbed and 
remains conformal along the flow, thus giving a control over the Liouville central charge. 
The partition function is evaluated numerically from combined analytic and perturbative 
information. In this paper we use the perturbative information up to third order. At special 
points the four-point integral can be evaluated and compared with our data. At the solvable 
point of minimal Liouville gravity we are in remarkably good agreement with the matrix 
model predictions. Possibilities to compare the result with random lattice simulations is 
discussed. 
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1. Introduction 

Liouville field theory (LFT) first appeared in the 2D quantum gravity context in the famous 
paper by Polyakov pQ. It has been demonstrated there that the LFT action appears naturally 
as the induced action of a conformal field theory living on a two-dimensional surface with 
non-trivilal background metric. Quantum LFT was recognized as a conformal field theory 
(CFT) and it is hardly an exaggeration to say that this peculiar example always stayed 
behind the motivations in the subsequent explosive development of the CFT, triggered by 
In particular, this induced LFT action governs a specific dynamics of the 2D quantum 
gravity, which is now called the critical (or Liouville) gravity. Characteristics of the critical 
gravity (critical exponents, correlation functions etc.) depend on the central charge of LFT, 
the last being simply related to the central charge of the inducing conformal matter. As a 
conformal theory LFT is rather complicated and differs in many important respects from the 
so called rational CFT's, developed and studied soon after the pioneering paper j2j. Despite 
essential progress in understanding of the space of states and hidden symmetries in LFT 
([HJ |3] or jlj), physical interpretation of the results in terms of fluctuating geometry of 2D 
surfaces remained unclear. 

Meanwhile, a very different approach to 2D quantum gravity has been under development. 
This approach can be called in general the discrete gravity. A smooth Riemannian manifold is 
substituted by a "random triangulation" , or more precisely, by a graph of the same topology 
(and possibly certain "matter" degrees of freedom attached to its vertices or faces) while the 
Riemann metric is simulated by the geodesic distances along the graph, the total number of 
vertices (or plaquettes) being interpreted as the surface area. The gravity problem is then 
reduced to a summation (in the spirit of the path integral) over the ensemble of such graphs. 
Progress in the study of these discrete models was drastically boosted after invention of a 
very powerful method to sum up the statistical ensembles of discrete gravity by reducing 
them to certain integrals over N x N matrices in the limit N — > oo [HI 13 El- The latter 
integrals can be sometimes evalutated explicitly, giving rise to exactly solvable models of 
discrete gravity. Despite the fact that not every discrete model can be treated by the matrix 
model approach, progress was so dramatic that now these discrete models are often called 
the matrix models of 2D gravity. In the limit N — > oo many types of critical behavior with 
different critical exponents had been discovered. Moreover, it is often possible to calculate 
explicitly the scaling functions which describe the "flows" between different critical points. 
Another impressive advantage of the discrete approach is that it allows to treat "surfaces" of 
arbitrary topology. This led even to the study of the grand statistical ensemble of surfaces 
with different topology, the object of extreme importance in the string theory. 

Approximately at the same time the Polyakov's gravity has been reconsidered in slightly 
different framework jU] (called sometimes the "Polyakov's" or "light-cone" gauge), which 
allowed to derive the critical exponents of this induced gravity in an explicit form (the KPZ 
scaling). It was soon recognized [TQJ [TT] that the same exponents are exactly reproduced 
in the original formulation of LFT provided certain finite renormalization of the coupling 
constant are conjectured and taken into account. 
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The nexus of the Liouville and discrete approaches to quantum gravity in fact flashed 
already in ref.jH] in the observation that some critical KPZ exponents are precisely those 
appearing in the critical behavior of random lattice systems. 

Since the very advance of the discrete 2D gravity [HI [Zj in 1985, it always challenged the 
field theoretical community. The ambition is not only to match with its achievements, i.e. 
learn to evaluate the correlation numbers and the scaling functions (and also to work at 
higher topologies). To reveal the field theory potential in 2D gravity to its full extent means 
to go beyond the limits accessible by the matrix models. In this paper I recall the most 
important steps and problems along this line. A scaling function for the spherical partition 
function of a perturbed CFT coupled to Liouville gravity is considered as a perturbative sum 
over the correlation numbers. Although it is still a problem to evaluate directly the whole 
set of these numbers, general arguments allow to establish important analytic properties of 
the scaling function. In order, these analytic properties restrict the scaling function to a 
large extent. In particular, sometimes even a few first orders of the perturbative series allow 
to reproduce the whole function, approximately but to a remarkably good accuracy. 

To my opinion the main lesson to draw from this simple exercise is that the Liouville 
gravity as formulated on the basis of early Polyakov's ideas is indeed a consistent and pithy 
structure whose physical content is rich and definitely extends beyond the relatively dull 
world offered by the matrix models in their present formulation ^2j (recent refreshments 
can be found in refs.jE]). 



2. Perturbed Liouville gravity 

Liouville gravity (LG) is the 2D gravity induced by a "critical matter". Its matter content 
is described by a certain "matter" CFT .Mm with the central charge cm- In the field theory 
approach LG is described as a joint CFT A4m+ Liouville + ghosts fj] 

£lg = £m + £l + £ g h (2.1) 
Conformal Liouville field theory is defined through the local Lagrangian 

£l = ^{d a <pf + /ie 2 ^ (2.2) 
Its central charge cl is related to the parameter b as 

c L = 1 + 6Q 2 (2.3) 

The combination 

Q = b- 1 + b (2.4) 

is called the background charge (basically for historical reasons). The scale parameter \i is 
called the cosmological constant because in LFT exp(2b(p)d 2 x is interpreted as the quantized 
area element of the 2D space. The ghost CFT is the usual BC system 

£ gh = -(CdB + CdB) (2.5) 

7T 
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of spin (2,-1) and central charge c g h = —26 [lj. The gravity interpretation requires the 
total central charge of the joint CFT ()2.1|) to be 0. In other words, parameter b in ()2.2|) is 
determined from the "central charge balance" equation 

c L + c M = 26 (2.6) 

Presently we are going to consider the perturbed Liouville gravity where the matter is 
not critical but described clS db perturbed CFT 

Anatter = + (2-7) 

Here $ is some relevant primary field from the spectrum of the matter CFT of dimension 
A and A is the corresponding coupling constant. It seems natural to take a simple solvable 
model of CFT as the unperturbed matter field theory, e.g., one of the famous minimal models 
Mp/ p i |2j. The perturbing field is then one of the relevant degenerate fields of this minimal 
model and its correlation functions are either known explicitly or constructible in principle. 
The matter central charge takes the worldwide known value c = c PiP / and therefore the 
Liouville parameters are rigidly fixed by ()2.6)1 . E.g., if the matter CFT is the free Majorana 
fermion CFT M.1/2 the Liouville parameter is a number b = a/3/4 and no more a parameter. 
For example, a question like "what is the physics of free fermion coupled to almost classic 
gravity?" cannot be addressed. From theoretical point of view it seems convenient to lift this 
restriction and get an independent control over the Liouville parameters. To this order let us 
add to the matter content of the theory another conformal matter of central charge c sp . This 
matter is unperturbed and always remain conformal. Therefore it formally does not interact 
neither with the Liouville mode (except for through the conformal anomaly) nor with the 
perturbed CFT (J2.7|) . Its main role is to contribute to the central charge balance ()2.6|) . For 
this reason we call it the "spectator" matter. To be precise, there is also some connection via 
conformal moduli. However in the simplest case of the spherical geometry (which we mainly 
address presently) any details of this "spectator" CFT except for its central charge are of 
no importance. It can equally well be considered as a formal shapeless ingredient seen only 
through c sp . In the last section we'll discuss possible realizations of the spectator matter. 
With this additional liberty the total matter central charge is 

c M = c + c sp (2.8) 

and the perturbed Liouville gravity is expressed symbolically through the Lagrangian 

jC x = £cft + £l + -Ue 2 ^ + £ sp + £ gh (2.9) 

As explained above, in the spherical case the last two terms in ()2.9j) are coupled only through 
the conformal anomaly. Therefore we will omit these background degrees of freedom, having 
in mind that they enter the central charge balance, which is solved for b as 



25 - c M _ /r-QM 

24 V 24 K J 
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For the branch chosen the semiclassical regime b — > is achieved when c sp — > — oo. In this 
limit the background space is in fact a classical symmetric sphere. Otherwise this parameter 
b characterizes the "rigidity" of the space with respect to quantum fluctuations. 

Interaction with the Liouville gravity "dresses" the perturbing primary field $ with an 
appropriate Liouville exponential field exp(2g0). Parameter g is determined from the fol- 
lowing kinematic requirement 

g(Q-g) + A = l (2.11) 
which means simply that the composite field 

U = $e 2g<p (2.12) 

is a (1,1) form and can be integrated over the space. Ordinarily the least of the two roots 
is chosen 

/ 25 - c M 1 1 -c M a , * 

9 = y—24 — y— +A (2 - 13) 

Notice that g is positive for relevant perturbations A < 1 (and negative for the irrelevant 
ones) . 



3. Spherical partition function 

The spherical partition function Z(fi, 0) of the Liouville gravity is known to scale as 

Z(/i,0) = /^/ 6 Z(l,0) (3.1) 

with some dimensionless constant Z(1,0). It is easy to see that the perturbing coupling 
constant A scales as p}l p with positive 

P = bg- 1 (3.2) 

The perturbed partition function Z(/i, A) introduces the scaling function Z of the dimen- 
sionless combination A/i _1 / p 

Z(fi, A) = /j Q/b Z {\^ l/p ) (3.3) 

This scaling function is the main object of our present study. Sometimes it will be 
more convenient to consider the "fixed area" partition function Za(X) defined as the inverse 
Laplace transform of Z(/i, A) 

Z A {\) = Aj^Z{^\)e» A (3.4) 

where the contour j goes along the imaginary axis to the right from all singularities of the 
integrand. In particular, in the critical gravity 
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Inversely 

ZM= / Z A (\)e-» A — (3.6) 

J(0) A 

where (0) at the lower limit indicates that the integral is often divergent at small A and 
()3.6|) means a properly regularized version. The name "fixed area" is motivated by the 
interpretation of the integral 

A— l e 2h<t> d 2 x 



as the total area of the gravitating sphere. For the same reasons as for Z(fi, A) the fixed area 
partition function has similar scaling form 

Z A {\) = Z A (0)z(h) (3.7) 

where it happens more convenient to normalize the scale invariant parameter h as 

/ a \ yp 

h = A (-J (3.8) 

4. UV expansion 

Partition function ()3.3)1 is a regular series in the coupling constant A 



Z(n,\) = Y t a n (-\) n (4.1) 



n=0 

Each term in this expansion is determined by the corresponding order in the interaction 
finH . The coefficients 

(27r)»n! J Mn K 1 y } 

are the integrated (and normalized) correlation functions of n insertions of the composite 
(j2.12j) . According to ()3.3|) the coefficients a n scale as a n ~ In the LFT approach the 

source correlation functions (i.e., the coordinate dependent functions before the integration) 
factorize into a product 

(U( Xl ) . . . U{x n )) = . . . $(x n )) M (e 2a ^ . . . e 2 ^«)) L (4.3) 

as it can be literally read off from the action (|2.9jl . The integral in (|4.2|) is not over the 
positions of all n insertions but rather over the n — 3 dimensional moduli space M n of the 
sphere with n punctures. Practically at n > 3 this amounts to place any 3 insertions of U at 
arbitrary fixed points xi, X2 and X3 and integrate over the positions of the remaining n — 3 
ones 

(U n )= f {W{x l )W{x 2 )W{x z )U{x i ) ...U{x n) )d 2 Xi ...d 2 x n (4.4) 
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$(x 1 )$(x 2 )) lA = \x 1 -x 2 \- 4A (4.6) 



Following the standard string theory practice we have also replaced the "nailed" insertions 
U(xx) by W(xi) = CCU(xi) to get rid of the redundant dependence on their positions and 
also to give the integrand the total ghost number 3 (as it is required by the ghost number 
conservation). Simplest example is the three-point function, where the moduli space is empty 
and the three-point function is simply the product 

(U 3 ) = (CC( Xl )CC(x 2 )CC(x 3 )) gh m Xl )<S>{x 2 )<S>{x 3 )) u (^V^Vs^s)^ (45) 

Two- and zero-point functions are simply figured out from the three-point one (see below). 

In the exactly solvable CFT models all multipoint correlation functions are in principle 
constructible. They are determined uniquely once the normalization of the primary fields is 
fixed. In what follows we use the standard CFT normalization of the matter fields through 
the two-point function 

($(xi)$(x 2 )) M 

The field normalization also gives precise meaning to the coupling constant A in (|2.7jl and 

(EH). 

In practice, even in solvable CFT's, it is usually difficult to go explicitly beyond the 
four-point functions (except for such special the free- field CFT's). At the same time 

the one-, two- and three-point functions are rather universal. Conventional CFT wisdom 
requires the one point function to vanish 

($(x)) M = (4.7) 

while the three-point function is fixed up to an overall constant C$$$ 

m Xl Mx 2 Mx 3 )) M = (4.8) 

1^12^23^321 

which is called the three-$ structure constant and is usually known explicitly in solvable 
CFT's. In (|4.8jl and often below we use the abbreviated notation xu = Xk — x\. Higher 
correlation functions have less universal form. The four-point case is more handy and allows 
explicit channel decompositions, e.g., 

m Xl Mx 2 Mx 3 Mx 4 )) u = Gm{X >% (4.9) 

F 14^23 1 

G M (x, x) = C'^C^iJ 7 ^ ^ ^ (Ailx^J 7 ^^ ^ |Aj|xJ 

in terms of the structure constants and the four-point conformal blocks related to the con- 
formal algebra with central charge c j2] 



x 7 (= 0) 
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These blocks are holomorphic in the projective invariant 

(x 1 - x 2 )(x 3 - ar 4 ) 



x 



(xi - x 4 )(x 3 - x 2 ) 



(4.11) 



The sum in (|4.9jl is over all primary fields Aj} entering the operator product expansion 
of$(aO$(0). 

The situation is nearly the same in LFT. The three-point functions are known explicitly 
for arbitrary exponential operators exp(2a</>) (the basic primary fields of LFT) 



<■ 



3 2ai$(a;i) 2a 2 4>(x2) 2a 3 <t>(x 3 ) 



Cx(ai,o 2 ,a 3 ) 



F12 



I2AX+2A2-2A3 



F23 



i2A 3 +2A 4 -2Ai 



X31 



2A 3 +2Ai-2A 2 



where 



A.; = a,i(Q- di) 



(4.12) 



(4.13) 



are the dimensions of the LFT exponential fields and the constant (letter a below stands for 
ai + a 2 + a 3 ) 



C L (ai, a 2 , a 3 ) 



n 



T fe (2a; 



T fe (a - Q) AJ- T fe (a - 2a. 



(4.14) 



is expressed in terms of certain special function T&(rr), which is an entire function related to 
the Barnes double gamma function. It's defining properties are the shift relations 



T b (x + b)= 7 (6x)6 1 - 2ba; T fe (x) 



(4.15) 



T b (x + b~ 



7 (6- 1 a;)6- 1+2 ^ lx T b (a;) 



It is also symmetric w.r.t. the reflection x —>■ Q — x and self-dual (i.e., invariant w.r.t. 
b — > 6 _1 ) as a function of the parameter b. Integral representation, which fixes also the 
overall normalization can be found e.g. in [T^]. In (|4.14jl and throughout the paper 7 (x) is 
the usual little 7 -function. The LFT four-point function 



^ e 2ai(f>(xx) e 2a 2 <j>(x 2 ) e 2a 3 <j>(x 3 ) e 2a,4(j>(x4)^ 

2(Ai+A 2 -A 3 -A 4 ) 



F24 



e 

2(Ai+A 3 -A 2 -A 4 ) 



(4.16) 



X34 



i4Ai 



F23 



2(Ai+A 2 +A 3 -A 4 ) 



0-2 0-4 



also allows holomorphic antiholomorphic decomposition [15] . similar to (]4.9|) 



x, a; 



ai a 3 

Q<2 G-4 

— C L ( 01, a 2 , — + ) C L ( — - iP, 03, a 4 ] J^p 



(4.17) 



Ai A 3 
A 2 A 4 



X jFr 



Ai A 3 
A 2 A 4 



x 
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In ()4.16|) and ()4.17|) Aj with % = 1, ... ,4 are the LFT dimensions (|4. 13)1 . The integral now 
is over a continuous spectrum of irreducible representations of the Virasoro algebra with 
central charge cl and dimension Q 2 /4 + P 2 . The corresponding four-point blocks 



Ai A 3 
A 2 A 4 



x 



Xj(=x) 



x 7 (= 0) 




x 3 (= 1) 



x 4 (= 




(4.18) 



are related to the conformal algebra with the Liouville central charge cl- 

Application of eq. (|4.4|) is not straightforward at n < 3. The moduli space has formally 
negative dimension and this is related to the non-trivial space of conformal Killing vector 
fields on the sphere with 2 and punctures. The simplest way to determine the correlation 
functions with n < 3 is to differentiate them sufficiently many times with respect to fi (each 
differentiation introduces an insertion of the area-element field exp(260)) until the order 
of the correlation function attains n > 3. This maneuver is quite a common way to treat 
the zero-mode problem of the sphere. In the case of the two-point function it's enough to 
differentiate once 







2! " , '- 28 *) Jj = C L ( g ,g,b)=(i r ry(P)) 
After an integration in \i 



(Q-2 9 )/b j(2gb - b 2 )~f(2gb 



-i 



TV fib 



<e 2 ^e 2 ^) L = - (7r/i 7 (6 2 )) 



(Q - 2s )/ 67 (2^-6 2 ) T (2^- 1 -&- 



tt(Q - 2g) 

This expression does not coincide with the Liouville reflection amplitude ?? 

< Q -2 g)/bl (2gb-b 2 ) 1 (2gb- 1 -b- 2 ) 



(4.19) 



(4.20) 



R(g) = -(n^(b 2 )Y 



(2g - Q)' 



(4.21) 



(which is ordinarily associated with the two-point function in LFT) differing in the factor of 
7r _1 (Q — 2g). Similar wisdom gives for the Liouville partition function 



<9/i 3 



C L (b,b,b) = -(7r^(b 2 )) 



Q/b 



(b 



-2 



7r 3 /i 3 &7(6 2 )7 (b 



and hence 



Q/b 



(4.22) 



(4.23) 



vr 3 Q 7 (6 2 )7(6- 2 ) 

Notice that this expression in not self-dual, unlike the three-point function itself. This fact 
looks somewhat surprising in the Liouville context and likely needs better understanding. 
Below we will meet many opportunities to support the consistency of ()4.20|) and ()4.23|) in the 
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2D gravity context. Therefore we have to take this breakdown of the Liouville self-duality 
in the 2D gravity application as the matter of fact. 

After all these preliminaries, we are in the position to construct explicitly few first terms 
in the perturbative expansion (j4.1j) 



a = l 

ai = (4.24) 
(^(fi 2 ))-^ Ql{2gb ~ ^^b- 1 - 6- 2 ) 7 (& 2 )7(^ 2 ) 



a 3 



8(2g-Q)(l-V) 
g, g) 



48vr 3 Z L 



9 9 



a4= 2M^I GMiX ^ )Gh \9 9 



x, x 1 d x 



In what follows we will be more interested in the fixed area partition function ()3.7|) . The 
scaling function z(h) has the following perturbative expansion 

z(/0 = 5>„(-/0 n (4-25) 

n=0 

where the dimensionless coefficients related to 

{ix^ n 'Pa n T{-l - b- 2 ) 



T{np- 1 - ft" 2 - 1) 

Explicitly up to the third order 



(4.26) 



*2 = (7(& 2 )) 

23 

In the last expression 



Zq = 1 

z x = (4.27) 
2g/b 7 (2^6 -6 2 ) 7 (6 2 )r(6- 2 ) 



8(6- 2 - 1)T(1 + 6- 2 - 2gb~ 1 ) 
C**$ Cl A=n \g,g,g) 



6(2tt) 3 z {A=n) 

L 



Cl {9 > 9 > 9) ~ T(3^ - Q)T^g)T(3gb^ - 1 - 6-2) (4 - 28) 
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is the fixed area A Liouville three-point function and the fixed area partition sum Z^ T reads 

z( A) _ f-_in\ Q/b r( 2 - b 2 ) 

L ~{ A J nW(b 2 )T(b- 2 ) [ ' 

For relevant perturbations ()3.8|) is small for small area A. This means that the pertur- 
bative expansion (|4.25|) in fact determines the behavior of the partition function at small 
surface area A. Notice also that the coefficients z n in ()4.26|) are suppressed at large n by 
additional T- function in the denominator as compared with a n . This means that if the 
perturbative series (j4.1j) is convergent, the fixed area z(h) is an entire function of h. 

5. Criticality 

From physical considerations we expect that at A —>■ oo the fixed area partition function 
behaves as 

Z A {\) = Z^A-Q'/t'e^ (5.1) 

where S vac (b) = —fi c is the specific free energy generated by the massive (finite correlation 
length) modes of the perturbed CFT Ai c - For dimensional reasons 

H Q = / (6)A" (5.2) 

with some dimensionless fo(b). This perturbed CFT interacts also with the quantized gravity 
and therefore S vac (b) reduces to the vacuum energy £ V ac(0) of the model (|2.7J) in "rigid" flat 
background only in the classical limit b 2 = 0. Otherwise it depends on b 2 and appears to 
be an important characteristic of the matter interacting with quantum gravity. Clearly the 
grand partition function Z(fx, A) develops a singularity at \i = /i c . Therefore S vac (b) is simply 
another interpretation of the critical point in the cosmological constant (usually referred to 
as the value of fi where the surface "blows up", one of the images convenient to make an 
illusion of understanding). The related singularity of Z(fi, A) in the coupling constant A at 



A -l#)J (5 ' 3) 

often happens to be the one closest to the origin. In this case it determines the finite 
convergence region of the perturbative series f)4.1j) . As it was already mentioned above, this 
in order implies that the function z(h) is entire. 

Let us also assume that the perturbed matter is "massive", i.e., the correlations of this 
coupled matter degrees of freedom do not extend far beyond the characteristic mass scale 
m c ~ A p . Then at A ^> m~ 2 the "massive" part of the matter "dies out", i.e., it does not 
contribute any more to the matter central charge. The gravity dynamics is determined then 
by the central charge of the "spectator matter" , which doesn't interact and therefore always 
remains conformal. Then the infrared parameters b' and 

Q> = l/fe' + b' (5.4) 
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are determined from the equation 



1 + 6(Q') 2 + c sp = 26 (5.5) 
These are the parameters entering the power-like part of the asymptotic (|5.1jl . More explicitly 



Q'= \iQ 2 + l (5-6) 



V 4 24 V 4 24 
In terms of h the asymptotic of the function reads 

log z(h) ~ 7r/o(6)A" + p (9 - 9L\ log/i + \og Zoo + . . . (5.7) 

where we have also added an undeterminate constant term. Dots stand for the subsequent 
decreasing terms in the large h expansion. In particular p is the order of the entire function 
z(h). 

Under certain additional assumptions about the analytic properties of z(h) the comple- 
mentary information provided by the asymptotic ()5.7|) and the first few terms in the series 
expansion ()4.25j) turns out quite restrictive and permits to achieve rather detailed description 
of this function. The efficiency of the combined analytic-numeric treatment implemented be- 
low depends decisively on these additional assumptions, which do not always hold true, being 
dependent on the choice of the matter field theory. From now on we restrict ourselves to 
a very particular example of the perturbed CFT, the scaling Yang-Lee model, where these 
additional properties are most likely valid. Therefore the central charge c in ()2.8j) and the 
perturbing dimension A in (j2.11|) are now fixed c = —22/5, A = —1/5. The "spectator" 
matter however remains undefined, the Liouville coupling b being a variable parameter of 
the model. 



6. Scaling Yang-Lee model 



The scaling Yang-Lee model 16 j is often used as a testing tool for various approximate 
methods in 2D field theory. It is probably the simplest (after the free field theories) example 
of perturbed CFT. The model arises as a perturbation of the non-unitary minimal model 
M.2/5- This CFT is rational and contains only two primary fields, the identity / of dimension 
and the basic field <p = $13 of dimension A = —1/5. The basic operator product expansion 
reads 

ip(xMO) = (xx) 2/5 [I] + {xxf'^C^ [^(0)] (6.1) 
where the three-y? structure constant is purely imaginary C wip = in with 

7 3 / 2 (l/5) 

/c = ' v ' ' = 1.91131269990474 . . . (6.2) 

57 i/ " : (3/5) 
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Below we will need also the matter four point function ()4.9|) 



(if(x)(f(0)(f(l)(f(oo)) YL = G Y l{x,x) = J r I (x)J r I (x) - k 2 J^(x)J^(x) (6.3) 

where the blocks are explicitly in terms of hypergeometric functions 

Fj(x) = x 2/5 (l - x) l/5 2 F 1 (2/5, 3/5, 6/5, x) (6.4) 
F v {x) = x 1/5 (l - x) l ^ 2 F 1 (l/5, 2/5, 4/5, x) 

The only relevant perturbation by the field ip gives rise to what is called the scaling 
Yang-Lee model. To make the theory real the corresponding coupling constant has to be 
taken pure imaginary. With this choice in ()2.9j) 

^matter = ^YL + ^V^)^ + £ S p (6.5) 

where g is determined by 1)2.11)1 . In a flat classical space-time this model is integrable. 
Its particle content and factorized scattering theory are known exactly The spectrum 
contains only one massive particle of mass m. Integrability provides us with a number of 
exact results. E.g., the exact relation 

m = k\ 5/12 (6.6) 
between the mass m and the coupling constant A is known 

5/24 



k ^ 



r(5/6)r(2/3) 



r(4/5)r(3/5) 
_ioor(i/5)r(2/5) 



1.2288903248... (6.7) 



The bulk vacuum energy is also obtained explicitly [121 E 



2 

m 



The classical space-time background is achieved in the limit 6 = 0. Hence the dimensionless 
characteristic fo{b 2 ) in ()5.2j) is known exactly at b 2 = 

/o(0) = ^= = 0.2179745... (6.9) 



7. Analytic-numeric treatment 

The analytic-numeric treatment here altogether follows the lines of ref.[2J) where precisely 
the same scaling Yang-Lee model has been considered on the rigid classical sphere. This 
rigid spherical background in our present case corresponds to the classical limit b 2 = 0. In 
this sense the present work is a direct extension of |2I] for the quantized geometry. 
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In [21] it was observed that the asymptotic ()5.7|) holds in the whole complex plane of 
h except for the negative real axis. All the zeros h n , n — 0, 1, ... of z(h) are real, all but 
the first one ho > being negative. In the present case of quantum geometry we make an 
assumption that these properties still hold for the whole allowed domain of the parameter 
< b 2 < b 2 . The upper limit here 

7 2 /6 

b\ = ~ = 0.4202041028 .. . (7.1) 



is related to the fact that at this value the "dressing" parameter g in (12.11)1 attains its critical 
value g = Q/2. For b 2 > b 2 the solutions to ()2.11|) are complex. The assumed location of 
the zeros is supported by the consistency of the numeric analysis of the next section. The 
presumption made, the asymptotic along the negative real axis reads (for brevity we suppress 
the argument b in the free energy f Q = fo(b)) 

z(h) ~ exp (71-/0 cos (vrp) (-h) p + . . .) sin (irf Q sin (717?) (-h) p + tt8 + ...) (7.2) 

where 

'-'(!-£)-§ 

The asymptotic location of zeros h n ctt n — ► 00 is hence estimated as 



/ sin(7rp) 

It is also convenient, as in to introduce the ^-function of these zeros 



^W=^ + ^p^ (7.5) 

The scaling function itself z(h) is related to (|7.5)l through 

log f^Mh^ (7.6) 

6 v ; 7 T sin(Trs) 2i7rs v ; 

where the contour goes along the imaginary axis to the left from the pole at s = 1 and to 
the right of the pole at s = p (note that p < 1 for the whole family of gravitational Yang-Lee 
models, see table El below). Inversely 

. sin7rs f 00 , .dlogzih — iO) 77 

Cyl (s) = / /i- s ^7 -<to 7.7 

where the contour of integration is shifted from the real axis to agree the branch chosen in 
(17.5)1 . By the way, from (|5.7j) and ()7.fij) it follows that 

resCvL(s) = fopsmnp (7.8) 



14 



and 



Cyl(O) =5 + 1/2 (7.9) 

Notice that 5 in (|7.3jl is always positive (see table EJ). Therefore ()7.3|) certainly doesn't 
give any approximation for the first zero ho, which, as it was mentioned above, is in any case 
positive. However for the next zeros, even for hi and h 2 , asymptotic estimate (|7.3|) gives 
surprisingly good approximations. This curious feature has been found out in the classical 
case b 2 = [21J and will be certified now for the quantum case by the numerical analysis as 
well as through the exactly solvable case b 2 = 0.4 (see below). This phenomenology motivates 
the following simple strategy, tried already in |21| : take ()7.4j) as a good first approximation 
for h n at n > and try to use the complementary information provided by the perturbative 
coefficients ()4.27j) to evaluate numerically the unknown parameter fo and also improve the 
approximation ()7.4|) as far as possible. 

To demonstrate this analytic- numeric procedure let's take the simplest version with only 
two first coefficients z\ and z 2 taken into account. The additional restriction for the zeros 
h n is contained in the sum rules 



1 

h Q ' ^— ' h r 



h 2 ^hl 



ra=l 
n=l 

Along the above strategy we approximate the sums in (j7.10J) through the asymptotic ()7.4|) 

_L_y ri = (7.11) 

1 + Y 2 r 2 = -2z 2 
K 

Up 



where Y = (/ sin(7rp)) IP and 



n = C(i/p, 1 - 5) (7.12) 
r 2 = C(2/p,l-5) 



Here ({s,a) is the usual Riemann zeta function. The system (|7.11j) is easily solved for 
ho = (Yn)- 1 and Y 2 = -2z 2 /(r 2 + r\) giving the approximation for f 



fo ] = —r^ ( ) (7-13) 

J0 sm(np) \r 2 + r 2 J y ! 

If first iV perturbative coefficients z\ . . . ,zn are known the corresponding iV sum rules 
()7.10|) can be solved for f and h n , n = 0, 1, . . . , iV — 2. The corresponding procedure, 
although quite simple, is described in appendix A. 



1 / -2z, 



p/2 
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Once the parameter Y and the first N — 1 zeros are restored in this order, the scaling 
function z(h) can be evaluated using (|7.6|) with the corresponding approximate Cyl(s) 



ins N ~ 2 1 

= Cyl(s) = B — + JZJ^ + F ^(^/P, N-l-S) (7.14) 

" n=l 



or 

N-2 



(7.15) 



8. Numerics 



In table C] numerical values of the non-trivial perturbative coefficients 22 and Z3 form ()4.27j) 
are summarized for different values of b 2 . At certain values of b 2 the four-point coefficient is 
also available and presented in the table. 



b 2 




z 2 


Z3 


z 4 


(sum) 


0.00 


00 


-0.0892857 


0.0229899 


-0.0318918 


-0.00318903 


0.01 


1.08 x 10 45 


-0.0883599 


0.0225977 




-0.00311130 


0.05 


7.43 x 10 8 


-0.0843801 


0.0209364 




-0.00278803 


0.10 


7909.85 


-0.0786500 


0.0186240 




-0.00235567 


0.15 


116.757 


-0.0718331 


0.0160092 




-0.00189425 


0.1974 


12.6331 


-0.0640680 


0.0132252 


-0.0014386452 


-0.00143865 


0.20 


11.4592 


-0.0635922 


0.0130616 




-0.00141311 


0.2404 


3.19037 


-0.0555991 


0.0104388 


-0.0010237712 


-0.00102381 


0.25 


2.47138 


-0.0534942 


0.00978831 




-0.000933453 


0.30 


0.799642 


-0.0409998 


0.00628732 




-0.000495034 


0.35 


0.327724 


-0.0255378 


0.00287061 




-0.000161857 


0.40 


0.155805 


-0.00714401 


0.000356752 


-8.50615 x 10~ 6 


-8.52261 x 10~ 6 



Table 1: Numerical values for the second, third, and sometimes forth order perturbative 
coefficients in the fixed area scaling function z(h). In the last coloumn we place the estimate 
of the four-point coefficient from the sum rules. 



Table |21 summarizes the results of the analytic-numeric procedure described in the previ- 
ous section. Here and are the estimates of the free energy parameter with the 
use of the perturbative coefficients from the table [T] up to z 2 , z 3 and also z A where the last is 
available. The exact values J^ exact ) come from the integrable classical case ()6.9)1 at b 2 = and 
from the matrix model solution in the case of minimal gravity b 2 = 0.4. These approximate 
values of the specific gravitational free energy fo{b 2 ) are plotted in figH Apparently it has a 
singularity at the "critical" value b 2 = b 2 . This singularity will be considered in more details 
elsewhere. 
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b 2 


P 


5 


A 2 ) 
Jo 


A3) 

Jo 


A \ 

Jo 


p (exact) 
JO 


0.00 


0.833333 


0.111111 


0.220407 


0.218156 


0.218036 


0.2179745 


0.01 


0.831646 


0.109935 


0.220318 


0.218091 






0.05 


0.824462 


0.106179 


0.219719 


0.217523 






0.10 


0.814333 


0.103698 


0.218260 


0.215950 






0.15 


0.802587 


0.103880 


0.215654 


0.213053 






0.1974 


0.789474 


0.106813 


0.211594 


0.208534 


0.208240 




0.20 


0.788675 


0.107063 


0.211308 


0.208218 






0.2404 


0.775255 


0.112194 


0.205864 


0.202209 


0.201858 




0.25 


0.771700 


0.113797 


0.204224 


0.200408 






0.30 


0.75 


0.125 


0.192522 


0.187677 






0.35 


0.719788 


0.142247 


0.171896 


0.165563 






0.40 


0.666667 


0.166667 


0.125625 


0.116905 


0.116151 


0.1158596 



Table 2: Specific energy f determined with the use of first two, three and four (when 
available) perturbative coefficients. 



To illustrate the convergence of the zero locations h n and also how a good approximation 
the asymptotic ()7.4|) gives even for h\, we collect some related numbers in table EJ The 
column corresponds to n — 1 in ()7.4|) with f = from table 121 The value of h^ xact ^ 
is taken at b 2 = 0.4 from the exact matrix model results (see next section) and at b 2 = as 
the best TCS fit of (21]. 



b 2 


h (2) 
a o 


a o 


a o 


, (exact) 




hf 




0.00 


2.42697 


2.43068 


2.43069 


2.43069 


-12.2410 


-11.7540 


-11.7667 


0.01 


2.44074 


2.44446 






-12.1991 


-11.7193 




0.05 


2.50261 


2.50648 






-12.0120 


-11.5656 




0.10 


2.60000 


2.60444 






-11.7880 


-11.3549 




0.15 


2.73096 


2.73652 






-11.5990 


-11.1519 




0.1974 


2.90514 


2.91258 


2.91264 




-11.4909 


-11.0075 


-11.0267 


0.20 


2.91685 


2.92443 






-11.4880 


-11.0020 




0.2404 


3.13562 


3.14598 


3.14606 




-11.5099 


-10.9730 


-10.9928 


0.25 


3.20130 


3.21256 






-11.5386 


-10.9867 




0.30 


3.69093 


3.70962 






-11.9506 


-11.2907 




0.35 


4.74483 


4.78283 






-13.3915 


-12.5134 




0.40 


9.22903 


9.38245 


9.38349 


9.38350 


-21.1991 


-19.3580 


-19.3879 



Table 3: Positions of first two zeros of h(z) at different approximations 



To summarize, using only three first coefficients of the power series in z(h) we have 
restored to an impressive accuracy the free energy / and also the first two zeros h and 
hi , the subsequent being given in this approximation by the asymptotic formula ()7.4|) . 
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With this data we can approximately calculate the next perturbative coefficients through 
the sum rules similar to (|7.10j) . Such estimate z^ nm ^ for the four point coefficient are quoted 
in the last column of table ^ These "extrapolated" numbers are compared with the actual 
ones where the last are available, i.e., at the classical point [21], the minimal gravity point 
b 2 = 0.4 and, in addition, at two special values of b 2 where the integral in the forth line 
of (j4.24[) is somewhat simpler and can be carried out (see in sect 11). The comparison 
is mainly to check the consistency of the four-point integrals with the expected analytic 
structure. Another reason is to estimate the numerical precision to be achieved in the four- 
point integration of (j4.24[) in order to add new information to the analytic-numeric data. 
The requirements to any numeric implementations turn out to be rather high (in most cases 
not less then 5-6 decimal digits of precision is needed). 

9. "Integrable" point b 2 = 2/5 

Point b 2 = 2/5 corresponds to c sp = 0, i.e., no spectator matter. The matter content of the 
critical (unperturbed) gravity is exhausted by a single minimal CFT model (M.2/5 i n the 
present case). This situation is called the "minimal gravity". There are many reasons to 
believe that the minimal gravity is solvable, at least in its "topological" aspects. Namely, 
all (or at least large classes of) integrated multipoint correlation numbers are in principle 
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calculable. The legend is that some versions of perturbed minimal gravity, to the utmost 
radical version all possible perturbations, are related to certain solvable matrix models 
(more precisely, the classes of their critical behavior). Anyhow, in many cases this folklore 
can be promoted to specific statements and our problem of the perturbed minimal Yang-Lee 
gravity is among them. 

The scaling function Z t (x,t), which describes the scaling region near the tricritical point 
of the one matrix model (at the planar limit corresponding to the spherical topology), is 
determined explicitly as 

d 2 

—Z t (x,t) = u(x,t) (9.1) 

through a solution u(x,t) of the following simple algebraic equation 

x = u 3 -tu (9.2) 

In ref. [22] this scaling region has been identified with the perturbed critical b 2 = 2/5 minimal 
gravity, parameters t and x playing the role of the cosmological constant and ^-perturbation 
coupling respectively. Of course, this identification doesn't fix the overall scale neither of the 
partition function ()9.1)1 itself nor of the coupling constants. Thus, we have to identify t and 
x with fi and A and the spherical partition function ()3.3j) Z t (x, t) up to certain normalization 
constants. 

Equations ()9.1|) and ()9.2|) result in the following expansion 

ZAx t) - t?' 2 V 1 r(3n/2 ~ 7/2) f — V (9 3) 

Z t {X,t) - t 2n , r(n/2 _ 1/2) ^3/2 ) 

„, , / 105 / x \ 105 / x \2 35 / x \3 

= z,(o, t )(i--( p75 ) + -( i57 j) --(^) +... 

The linear in x term is regular and therefore cannot be described unambiguously in the 
field theoretic context. In particular, this term doesn't contribute to the fixed area partition 
function below. Hence we do not pay much attention to this term and instead compare the 
next terms with the perturbative coefficients ()4.24j) . The x 2 term leads to the relation 

x_ _ Xl eg 

£3/2 (^3/2 ^" > 

which is also consistent with the x 3 term. Here 

l es = = 0.0845223 . . . (9.5) 

Also the singular point in t is determined from (|9.2j) and (|9.4j) 

3 (I \ 2/3 

^ = -[y) = 0.11585962159187 (9.6) 
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This is the number quoted as [ n the table El 

Neglecting the overall normalization of the partition function, which in any case doesn't 
enter the fixed area scaling function z(h), we are in the position to write the latter explicitly. 
The coefficients in (I4.25J1 are 

2v^F(-/ cg )" 



Numerically 



n\Y{n/2 - 1/2) 



(9.7) 



Z-2 



z 3 



2,1 



-0.00714401 



0.000356752 



(9. 



6 eg 



-8.50615... x 10" 



The last number has been used in the previous section for the numerical analysis. 
The following integral representation 

z(h) = —j= ^y 2 dyexp (y 2 + (9.9) 

where the integration is parallel to the imaginary axis and goes to the right from the singu- 
larity at y — 0. 

In figEl exact scaling function ()9.9)1 is plotted. For comparison we present, without any 
comments, the numbers obtained by the order 3 (i.e., first three coefficients zi, z 2 and z 3 are 
used in the analysis) analytic-numeric procedure of sect. 8. 

It is worth mentioning that the family of gravitational Yang-Lee models < b 2 < b 2 
contains another interesting point b 2 = 0.3, where it can be reduced to the minimal grav- 
ity and therefore there are all reasons to expect the model to be solvable. At this point 
c sp = —22/5 and the model can be thought of as two independent Yang-Lee models, one is 
perturbed while the other remains conformal (spectator). So far nothing special, one can 
make many similar combinations of minimal matter models constructing different theories of 
gravity, most of them being non-solvable (or non-integrable, in a wide sense). Combination 
•M.2/5 © M- 2/5 is exceptional, because it admits an alternative description in terms of single 
minimal model M.3/10 (central charge —44/5) [2H]. The table of degenerate dimensions of 
this model 



2 


49/40 


3/5 


1/8 


-1/5 


-3/8 


-2/5 


-11/40 








-11/40 


-2/5 


-3/8 


-1/5 


1/8 


3/5 


49/40 


2 



(9.10) 



includes the fields $1,5 and $2,5 of dimension —1/5. In the non-diagonal (.D-series) version 
of this minimal model certain combination of these operators is interpreted as ip acting only 
inside one of the constituent M.2/5 models. It is natural to expect that such perturbed 
minimal gravity can be related to an appropriate matrix model, and therefore is solvable. 
The point seems to merit further study. 
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Figure 2: Continuous line: exact function z(h) as defined by Open circles: the 

canonical product over the zeros recovered from Z2 and z^ through the analytic-numeric 
procedure 



10. Four-point integral 

Relative success of the numerical program of sect. 8 makes it tempting to improve the approx- 
imation by involving the next four-order term z±. This coefficient is related to the four-point 
correlation number. In the framework of Liouville gravity it requires an integration over 
moduli f)4.4j) of the product of the matter and the Liouville four-point functions, as in the 
last line of eq. (|4.24jl . In the fixed area context it is more natural to use the fixed area 
Liouville four-point function 

n (A=ix)f _x G L (x, x)\ x rmi\ 

°- {x > x) = T(4 gb -i-Qb-i) (101) 



The fixed area term is rewritten as 



z A = (a ^ 10 - 2 

24(2n)%Z { L A=n) 
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where, as usual 



and 



, {A=n) 7 Q/ fc (& 2 )r(2 - b 2 ) 

' L nW(b 2 )T(b- 2 ) 



1(b) =b J G YL (x,x)G 



7T \x,x)d 2 x 



(10.3) 
(10.4) 



with Gyl(£,^) from eq. (|6.3j) . 

In general the Liouville four-point function admits the block representation (|4.17|) . which 
reads in our symmetric 



G^x,x)\^ =1 = n g f^r g (P)f P ( * 



9 A 9 



Here the overall factor lZ g reads 



while the function 



K g =h(b*)b 



,(P) 



-n 2 T\(2g - Q/2) 
7r 2 T 6 (2zP)T 6 (-2zP)T£(2 5 - Q/2) 



(10.5) 



(10.6) 



^ ' Tl(b)Tl(2g - Q/2 - iP)T 2 (2<? - Q/2 + zP)T 2 (Q/2 - zP)T 2 (Q/2 + iP) 
admits the following integral representation (convergent at g > Q/4) 

dt sin 2 Pt (cosh 2 (Q - 2g)t - e~ Qt cos 2 Pt) 



10.7) 



r 9 (P) = sinh27rfe" 1 Psinh27rfePexp -8 



t 



sinh 6t sinh b _1 t 



(10-8) 

Integral f)lU.5j) should be taken literally if Q/2 > g > Q/4, otherwise there are additional 
"discrete" terms due to the singularities of the integrand The prime near the integral 
sign in (jlO.Bj) indicates this subtlety. 

Apparently the form Gyl{x, x)G^~ n \x, x)d 2 x in the integral ()10.4j) is invariant under 
the group of 6 modular transformations 

R : x — > 1 — x 
I : x — > 1/x 



RI :x^ 1/(1 -x) 
IR : x — > 1 — 1/x 
IRI = RIR : x -> x/(z - 1) 



(10.9) 



Therefore the integration in (|10.4|) can be restricted to its fundamental domain F ={Rex < 
1/2; |1 -x\ < 1} (see figEJ) 



1(b) = 6b [ G YL (x,x)G[ A - 7T \x,x)d 2 x 



10.10) 
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Evaluation of integrals in ()10.4j) or (jlO.lOj) . even numerical, presents certain technical 
difficulties. They are mostly related to the complicated nature of the Liouville four point 
function. Expression ()10.5p involves additional integration in (jlO.Sj) to evaluate the special 
functions. This makes the integral in (jlO.lOj) three-fold. In addition the general conformal 
block entering (jl0.5|) is usually known only as a power series in x, its evaluation requiring 
certain recursive procedure [21]. 

There is another representation of integral (jl().4j) which takes advantage of the holomor- 
phic factorization in each term of the matter and Liouville block decompositions ()6.3|) and 
(j4.17|) respectively. It is based on the following simple calculation. Consider the integral 




(10.11) 



where G(x, x) is singular at x — 0, 1 and oo only and admits the representation 

G(x, x) = J2 AiFi{x)Fi(x) (10.12) 
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(the sum over % might be as well an integral over a continuous spectrum or both) with some 
coefficients A4. The holomorphic "blocks" are singular at x — 0, where J-i(x) ~ x Ai , and also 
at x = 1 and x = 00, the combination ()10.12j) being arranged so as to form a single- valued 
function of x. Straightforward application of the Stokes theorem gives 



1 - Yi^ AiIiJi 



(10.13) 



where 



J; 



Ti(x)dx 
Ti{x)dx 



(10.14) 



c 



and the contour C goes from —00 to —00 around x = counterclockwise. If, in addition the 
integrand in (|10.11|) is symmetric with respect to the modular group ()10.9|) . equation ()10.13J) 
is reduced to 

I = - ^ shM> A i) Ail? (10.15) 

i 

In our particular problem this identity reads (again up to certain discrete terms at g < 
0/4) 



1(b) 



-bn n 



I" dP 

f I^( P ) (sin7rA P /J(P) - K 2 sin7r(A P - l/5)iJ(P)) (10.16) 



where 



h{P) 



T lf (x)J r P 



A, 


A, 


A, 


A, 


A, 


A 3 


A 9 


A 9 



a; I dx 



10.17) 



x I dx 



J-i(x) and J-tp(x) are from eq. (|t>.4j) . 

Expression f)lU.16j) has an obvious advantage compared with the form (jlO.lOj) . since it 
involves one integration less. However, evaluation of the line integrals ()10.17j) requires the 
Liouville block to be computed to the required accuracy up to the upper limit x — 1. For 
a function known mainly through its power series this makes an additional problem. On 
the other hand, representation ()10.10|) contains only the blocks with \x\ < V3/2, where the 
convergence of the block power series is still reasonable. To conclude, it is difficult to say 
in advance which one of these two representations is better for numerical work. In the next 
section we take two special values of b 2 where the Liouville correlation function can be found 
in closed form. The Liouville blocks explicitly known make formula (jlO.lfij) very effective 
and accurate. This will allow us to estimate the numerical precision of the direct integration 

in pm 
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11. Special b 2 

Liouville correlation function ^e^i^i) . . . e 2a N <f>{x N )^ Qi eve i p poles if its parameters a« satisfy 
certain relations. Following A.Polyakov j2S] we will call them the "on-mass-shell" conditions 
and the corresponding correlations the "resonant" ones. The simplest on-mass-shell condition 
reads for a = Yli=i a i 

a = Q-nb (11.1) 

with n non-negative integer. It can be easily traced to the divergency in the integration over 
the "zero mode" of the Liouville field (2SII2ZI- The residues in the resonant poles 1)11.1)1 can 
be "naively" red off from the Liouville Lagrangian (J2.2)) 

( p 2 ai <t>{xi) 2a N <t>{x N )\ . , (-/i) W gn(«l,---,ajv) (WA 

\ • n! a-Q + nb 1 } 



where 

Gnicti, ■ ■ ■ , a 



n ) = j ^ e 2ai#Cl) . . . e 2a ^ iXn) f[ e 2H{yk) d 2 y^J (11.3) 

and the expectation value in the right hand side is over free massless boson without the zero 
mode. 

In particular in the four-point case the simplest n = resonance doesn't require any 
integrals in the right hand side and the "reduced" four point function of (J4.16)) takes the 
form 

a 2 a 4 I a - Q 



Gi 
where 



Go(x,x) = (xx)- 2aia2 ((l - x)(l - x)Y Iaias (11.5) 
Similarly, the next resonance at a = Q — b is controlled by a single "screening integral" in 

dm 



Qx(x,x) 



1 f #y 

(xx) 2aia2 ((l -x)(l - x)) 2a ^ J (yy) 2a > b ((l -y){l- y)) 2a * b ((x - y){x - y)) 2a ^ b 



11.6) 

It admits the two-term holomorphic decomposition 

Qi(x, x) = kiFi^F^x) + k 2 F 2 (x)F 2 (x) (H-?) 

where 

k _ 7r 7 (l-2o 2 6-2oi6)7(l-2a 3 &) fn g) 

1 7(2 - 2a 2 b - 2a 3 b - 2 ai b) [ ' ' 
= 7r 7 (l-2a 2 &)7(l -2 ai b) 

2 ~ 7 (2 - 2a ± b - 2a 2 b) 
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and 

F x (x) = x- 2aia2 (l-x)~ 2aia3 2 F 1 (2a 1 b + 2a 2 b + 2a 3 b- 1, 2 ai b, 2a x b + 2a 2 b, x) (11.9) 

F 2 ( X ) = x l-2a 2 6-2a 1 6-2a 1 a 2 ^ _ x y^a z _ 3^ 2 a 3 6, 2 - 2di& - 2a 2 6, s) 

To check the consistency of the resonance condition in LFT, in Appendix B we rederive the 
residues (111.5)1 and (jll.7j) from the general four-point expression (|4.17j) . 

From eq. (|l().l|) it is clear that the poles at pi. 1)1 are canceled out in the fixed area 
correlation functions and the finite parts 

( e 2ox*(*i) _ _ _ e 2a N ^x N )^(A) = b -l A -ng n ( au . . . ? (11. 1Q) 

are expressed through pi.3j) . In our four-point case this means that p().4)l reads 

Z(b n ) = — I G YL (x,x)Q n (x,x)d 2 x (11-11) 



7T 

1. The first resonant point n = appears in our particular problem pi. 4)1 at b = bo 

H — 4,/f; 

b 2 = — = 0.240408205773 . . . (11.12) 

At this point the Liouville dressing parameter g 2 = 2/5 so that 

g (x, x) = [xx{l - x)(l - x)]~ A/5 (11.13) 

Hence 

[xx(l — x)(l — x)] 4//5 

Direct integration over the fundamental region gives (NIntegrate of Mathematica 3.0) 



J(6o) = 6 / „ GYL l*' X) ^ /5 d 2 x = -59.9029678... (11.15) 
Jf \xx(l — x 1 — x)\ 



'f [xx{1 — x){1 — x)\ 
The reduction formula p0.15)) gives in this case two terms 



I(b ) = sm—I 2 -K 2 sm—I 2 (11.16) 



Here both contour integrals are carried out explicitly 



/■ 1 F(2/5, 3/5, 6/5, t) _ 2^n(V5 - l)r 2 (6/5) 

1 J *2/5(i_t)3/5 r(7/io)r(9/io)r(9/5) 1 ' ] 



^(1/5,2/5,4/5,*)^ _ ^(1/5) 



o t 3 /5(i _ t)3/5 2 2 / 5 r(7/10) 
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and ()11.16|) sums up to 



J(6o) = _^- 7 5 (i/5) = -59.9030606... (11.18) 



Finally 

z A = — J ^ ... = -0.00102377123. .. (11.19) 

3846tt 4 ^ A) 

is the number used in sect. 8 for this special point. 

2. The second resonant point n — 1 corresponds to 

6?= — = 0.197368421... (11.20) 
76 

2 19 

q = — 

y 60 

Notice that gb = 1/4 and thus the integral (jll.6|) reads 

1 f d 2 y 

(x ' x) = [(xx)(i-x)(i-x)p a y [^(i-,)(i-^)(*-,) ( *-^ (11 - 21) 

This fact is general for the n — 1 resonance in the symmetric case, for Ylt=i = Ag = Q — b. 
The integral is a degenerate case of (jll.6|) . The result is recovered at the limit 2oj6 — > 1/2 
in (fTTTjl 

c , -x 4[K(x)K(l-x) + K(l-x)K(x)] 

^ = [XX(1-X)(1-X)]W30 ( 1L22 ) 

in terms of the elliptic integral of the first kind 

K(x) = - C . dV = (11.23) 
2 Jo y/y(l-y)(l-xy) 

Collecting all together we find 

Ifa) = ^ J G YL {x,x)g 1 (x,x)d 2 x (11.24) 

Direct integration gives with Mathematica 3.0 

6 f 

X(bi) = - Gy L (x, x)g x {x, x)d 2 x = -302.01932532722 . . . (11.25) 

The reduction formula ()10.15j) doesn't apply here directly due to the degeneration. It is more 
convenient to reduce first a more general integral with ()11.6|) as the Liouville part. Then the 
limit 2aj6 — ► 1/2 reveals 

q / pj X37T "\ 

1(h) = - I sin ^-IjJj - 2k 2 sin — (11.26) 

7T \ dU OU / 
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where 

In numbers 
and 



'I* l Q _ s ))19/30 C!a: > ^ J Q _ x)) 19/30 dX 



X(bi) = -302.01932904124 . . . (11.28) 



z A = — J ( &1 ) = -0.0014386452412875 . . . (11.29) 

3846vr 4 Z{ A) 

Reduction formula gives more accurate numbers. This permits to estimate the precision 
of the direct integration. In the above two examples it varies between 5-8 decimal digits. The 
accuracy of the direct integration is determined mainly by the singular behavior near x = 
in the fundamental region. Thus, for b 2 closer to the critical point b 2 , where the "tachyon" 
divergence appears, the direct integration is hardly expected to give a good approximation. 

The following remark is in order here. At the fixed area picture the n-th resonance 
(jll.lj) correlation function gives in the partition function a contribution proportional to the 
negative integer powers A~ n . The transform to the grand partition function (|3.fij) results in 
the terms containing fj, n log \x. This is quite usual phenomenon in the field theory: integer 
powers of the coupling constant are usually decorated by logarithm. Otherwise these terms 
are regular in /i, i.e. contact (from the common point of view they are out of the field theory 
scope). On the other hand, the analysis of the critical behavior in solvable matrix models is 
similar to that in the mean field, or Ginzburg-Landau theory, and thus deals only with the 
algebraic functions where there is no place for the logarithms. This might give an impression 
that the mean field like pattern of critical behavior is common for the statistical systems 
interacting with quantum gravity. We have seen in the above examples that this is certainly 
not the case and might be true at most for the systems related to the solvable matrix models. 



12. Discussion 

I think that the main lesson to be learned from those rather random calculations presented 
above is that the perturbed Liouville gravity based on the existing LFT constructions is 
indeed a consistent tool to treat the standard problems in 2D gravity. Although at present 
it cannot yet give any exact description even in the problems known to be exactly solvable, 
sometimes it works equally well in the situations which are not covered by the solvable 
matrix models and therefore most probably are not exactly solvable. On the other hand, 
non-solvable gravities certainly exist and might be observed e.g., as critical points of non- 
solvable random lattice statistical systems. More generally, it is likely too hasty to claim the 
field theory obsolete, substituted by the matrix models (or whatsoever). 

As for the "spectator matter" appeared in an apparently artificial way, there are many 
practical ways to introduce such content in the statistical models of lattice gravity. The 
simplest way is to add a Z)-dimensional free massless lattice boson (without zero mode) 
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which doesn't interact with other matter degrees of freedom and contributes simply as an 
extra weight factor det~ D ^ 2 A lat , where A lat is the discrete Laplace operator of the graph. 
Apparently this parameter D = c sp can be as well made continuous. Not to talk about the 
obvious idea of putting several extra non-interacting spin-like degrees of freedom and tuning 
them to their critical points. As far as I understand, all these situations can be studied either 
through the extrapolation of the finite lattice ensembles or in the Monte-Carlo simulations. 
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A. Series analysis 

Suppose that the order p of z{h) is less then 1 and therefore 

n=0 ' n=0 \ k=l / 

where A" is any natural number (for our particular problem N > 2) and 



* fat 



Thus 

P N (h) = z{h) exp (f^rki-hA (A.3) 

is a polynomial of order N — 1 with zeros at h n , n = 0,1, . . . , N — 2 (we have introduced 
certain unnecessary minus signs in the notations for further convenience). 

This set of trivial identities gives rise to the following algorithm for the problem of sect. 7 
: solve first truncated A" sum rules (17.10)1 for Y and first A^ — 1 zeros ho, ... , given the 
first A" non-trivial coefficients z„ in 



:(h)=J2zn(-hT (A.4) 



n=0 
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The truncation means that all terms in the r.h.s. of (|7.10|) with n > N — 2 are replaced by 
the corresponding asymptotic values (|7.4J1 . In our present algebraic context the coefficients 
r k in (|A.2J) are taken as 

r k = d k Y k (A.5) 
with some known numbers d k (in our actual context d k = k~ l C > {k/ p, N — 1 — 5)). Then 

(oo \ oo 

^^HM =Y^a n (-hY) n (A.6) 
fe=l / n=0 

is a regular series in hY with known coefficients 

«■= E n| < A -7) 

fcl+...+fe;=n £ 

The coefficient of order A" in (|A.3|) is zero. This gives the algebraic equation for Y 

N 

Q N (Y) = Y,akZ N - k Y k = (A.8) 

fc=0 

Once Y is found as a suitable solution to this equation, the A^ — 1 roots of the polynomial 
are found as the roots of the polynomial 

N-l 

52Q n (Y)(-h) n = (A.9) 

n=0 

where 

n 

Qn(Y)=J2 a kZn-kY k (A.10) 
k=0 

B. Resonant Liouville correlations 

Consider the general Liouville four-point function (j4.17|) 

A c ;x J T 



Ai A 3 
A 2 A 4 



A,.; x 



in the vicinity of the first resonant point 

o-i + a 2 + a 3 + a 4 - Q = e (B.l) 
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with e small. Due to the denominators T&(ai + a 2 — c) and T fe (a 3 + a 4 + c — Q) in 



C L (oi, a 2 , c)C L (Q - c, a 3 , a 4 ) : 

T 6 (6)T ;) (2a 1 )T 6 (2a 2 )T ;) (2c 



/ o> o of,2\ (Q-ai-a 2 -a 3 -a 4 )/b 



X 



T 6 (oi + a 2 + c - Q)T 6 (ai + a 2 - c)T 6 (c + a x - a 2 )T fe (c + a 2 - 
T 6 (fc)T 6 (2a 3 )T 6 (2a 4 )T 6 (2Q - 2c) 



(B-2) 



T 6 (a 3 + a 4 + c — Q)T b (a 3 + a 4 — c)T fe (c + a 3 — a 4 )T b (c + a 4 — a 3 ) 

there are close poles: at c = ai + a 2 — e to the left from the integration contour and at 
c = ai + a 2 to the right. This pinch results in the singular contribution 



2e 



A x A 3 
A 2 A 4 



A 



«l+«2 ! 



A x A 3 
A 2 A 4 



A ai 



+0.2 ) 



(B.3) 



All other multipliers in ()B.2|) cancel out in the residue. It is not a problem to verify that 



T 



Ai A 3 
A 2 A 4 



A. 



„— 2aia2 i 



l-x) 



(B.4) 



A "mirror" contribution comes from the pinch at c = Q — a\ — a 2 provided by the T;,(ai + a 2 + 
c — Q) and T^(a 3 + a 4 — c) denominators in (jB.2|) . Its contribution is identical and we arrive 
at (|11.5j) . Of course the same situation appears if one or more of ai in (|B.1J) are reflected, 
e.g., a\ + a 2 + a 3 — a 4 = e or Q + ai — a 2 + a 3 — a 4 = e. Then other pairs of T;,-functions 
give rise to the singularity, which remains the same up to the expected reflection factors. 
Let's now take the second resonant point 



a% + a 2 + a 3 + a 4 - Q = -b + e 



(B.5) 



The same pair of denominators T 6 (ai + a 2 — c) and T 6 (a 3 + a 4 + c — Q) in (jB.2|) produces 
two pairs of close poles at Ci = a x + a 2 , a x + a 2 — e and c 2 = a x + a 2 + 6, a x + a 2 + 6 — e. 
The singular residues are 



7r/i 7(1 — 2ai& — 2a 2 6) 
"2e 7(2a 3 6)7(2a 4 6) 



J 7 



Ai A 3 
A 2 A 4 



A 



«l+«2 1 ' 



and 



7r/x 7(1 — 2a 3 6 — 2a 3 b) 
27 7(2a 2 6)7(2a!6) 



J 7 



Ai A 3 
A 2 A 4 



A 



ai+a 2 +6i 



Ai A 3 
A 2 A 4 



Ai A 3 
A 2 A 4 



A 



«l+«2 1 



A ai + a2 +6, x 



(B.6) 



(B.7) 



The "mirror" pair T&(ai + a 2 + c — Q) and Tb(a 3 + a 4 — c) contributes identically (in fact, 
this is a general feature related to the symmetries of the integrand in (j4.17j) ). The identities 
for the blocks 



x 



Ai A 3 
A 2 A 4 



A ai + Q2 , X 



fB.8) 



x 



- 2oia3 2 Fi(2ai6 + 2a 2 6 + 2a 3 6 - 1, 2a x b, 2a x b + 2a 2 b, x) 
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and 



T 



A ai -(-a2+6i 



A 2 A 4 

x l-2a 2 6-2a 1 6-2a 1 a 2 ^ _ ^'^a^p^ _ 3^ 3^ 3 _ 3^ _ 2 a 2 6, X ) 



(B.9) 



bring us to eq. (|11.7jl . 
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